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ABSTRACT 

We revisit the properties of white dwarfs accreting hydrogen-rich matter by constructing steady-state models, 
in which hydrogen shell burning consumes hydrogen at the same rate as the white dwarf accretes it. We obtain 
such steady-state models for various accretion rates and white dwarf masses. We confirm that these steady 
models are thermally stable only when the accretion rate is higher than ~ IQt'Mq yr"'. We show that recent 
models of "quiescent burning" in the "surface hydrogen burning" at a much wider range of accretion rates 
results from the too large zone mass in the outer part of the models; hydrogen burning must occur in a much 
thinner layer A comparison of the positions on the HR diagram suggests that most of the luminous supersoft X- 
ray sources are white dwarfs accreting matter at rates high enough that the hydrogen burning shell is thermally 
stable. Implications on the progenitors of Type la supernovae are discussed. 

Subject headings: accretion - binaries: close - stars: evolution - novae, cataclysmic variables — supernovae: 
general - white dwarfs 



1. INTRODUCTION 

It is widely accepted that progenitors of Type la super- 
novae (SNe la) are mass accreting C-O white dwarfs, although 
the nature of the progenitor binary systems is still under de- 
bate (e.g.,[No moto et al. 1997; Hillebrandt & Niemeyer 20001 
iNomoto et alJl2003l) . In the most popular model, the Chan- 
drasekhar mass model, a C-O white dwarf accretes mass until 
it gro ws in mass to M = 1 .38 Mf?^, then explodes a s an SN la 
(e.g., Nom otoll982HNomoto et al]|T984l;IUvioll200 0). To find 
a path to SNe la is to find a way for a C-O white dwarf to 
increase its central density and temperature to carbon ignition 
as the Chandrasekhar limit is approached. If the mass donor 
is a normal star, the stability of the hydrogen burning shell in 
the accreting white dwarf is crucial for its evolution. 

A hydrogen-shell burning is ignited in an accreting 
white dwarf when a certain amount of hydrogen-rich mat- 
ter is accumulated in the envelope. The shell burning 
is unstable to flas h if the accretion rate is lower than a 
critical rate (e.g., ISugimoto & Mivaiil Il981h . The flash 
is stronger if the white dwarf is coole r and the accre- 
tion is sl ower (iPaczynski & Zyflii:owlll978t ISion et al.1 [19791: 
lFuiimotolfT982allbh . If the shell flash is strong enough to 
trigger a nova outburst, most part of the envelope would 
be lost from the system. Moreover, a part of the orig- 
inal white dwarf matter could be dredged up and lost in 
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the outburst wind. Therefore, the white dwarf mass could 
decrease aftei^ one cycle of nova outburst (e.g., Prialni^ 
1986; Prialnik & Kovetz 1995; lYaron et alj I2005i) . although 
iTownsley & Bildsten C2004.) did not find such indication in 
the nova systems. 

If the accretion rate M is high enough, on the other hand, 
the shell burning is stable and the mass of the white dwarf in- 
creases with accretion. Thus, the critical accretion rate Mjtabie, 
above which hydrogen shell burning is stable, is an impor- 
tant physical value for the evolution of accreting white dwarfs. 
(For M being appreciably higher than MstaWe, see §5). 

Studies on the t hermal instability of a th i n nucl ear-burning 
shell began when ISchwarzschild & HarmI (Il965h found that 
helium shell flashes were caused by a thermal instability of 
a shell around a C-O core of an evolved star. Since then, 
the stability of a nuclear shell burning in an accreting white 
dwarf has been extensively studied in connection with nova 
outbursts ("e.g.. ISugirnoto & Fuiimotolll978l and references 
therein). ISienkiewiczl (Il980h obtained the stability bound- 
ary by examining thermal stability of steady-state models 
for accreting white dwarfs of various masses. The stabil- 
ity boundary, which is MstaWe — lO"^M0yr"' at M = O.8M0 
and higher for more massive white dwa rfs, is consistent with 
fully time-dependent calculations (e.g., Paczvhski & Zvtkow| 
1978; Sionetal. 1979; Iben 1982; Prialnik & Kovetz 19951 
l^ron et al.ll2005l). 

Recently, IStarrfield et all (l2004l) published new models of 
accreting white dwarfs called "surface hydrogen burning" 
models. They tried to show that if a very hot white dwarf 
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after a nova outburst started accreting hydrogen-rich matter, it 
developed surface hydrogen burning that stably converts hy- 
drogen into helium and helium into heavier elements; then 
the white dwarf mass could grow to the Chandrasekhar limit. 
They obtained such "surface hydrogen burning" models for 
accretion rates ranging from 1 .6 x 10"^ to 8.0 x lO'^M© yr"' , 
and identified these models as supersoft X-ray sources. The 
properties of their models, however, clearly contradict previ- 
ous results of evolutionary and steady-state model s (e.g.. llbenl 
[1981 iPrialnik & Ko"vet3[l99l ISienkiewiczllT 98g). 

We will show in this paper that the contrasting proper- 
ties of the "surface hydrogen burning" models result from 
too coarse zoning for the envelope structure. We first re- 
calculate steady-state models for various accretion rates and 
white dwarf masses. If accretion continues long enough, 
the mean luminosity of the white dwarf should approach the 
steady-state luminosity irrespective of the initial luminosity. 
We then examine the thermal s tability of these models to con- 
firm the results of lSienkiewiczi ( 1980) using the updated input 
physics. 

In obtaining our steady models, we solve the full stellar 
structure in o rder to avoid the a pproximations made in the one 
zone models (lFujimotoll982ah . We present more detailed nu- 
merical models, their stabilities, and the exact stability bound- 
ary values in tabular forms and formula for a wider range of 
the white dwarf m ass and the accretion rate compared with 
LSienkiewiczl ([T97i ri980). 

We have adopted linear stability analysis mainly because 
we can avoid long evolution calculations for many mode ls. In 
Oder to judge the stability by evolution calculations as in llbeiil 
^982), we have to adopt sufficiently short time steps, which 
demands long cpu times. Also, by a linear analysis we can 
determine the stability boundary exactly, while it is difficult 
for the evolutionary calculations, because near the stability 
boundary the growth time is so long that it would be difficult 
to distinguish the variation due to thermal instability from a 
normal evolutionary change. 

Then we examine the "surface hydrogen burning" models 
by IStarrfield et al.l (l2004t) and discuss the reason for the dis- 
crepancy with the previous results. Our numerical methods 
are described in §2 and our results are given in §3. In §4 we 
discuss the reason for the discrepancy between Starrfield et 
al.'s (2004) models and the previous results, and compare the 
loci on the HR diagram of our steady-state models to those of 
luminous supersoft X-ray sources. Implications on the pro- 
genitors of Type la supernovae are discussed in §5. 

2. STEADY-STATE MODELS 

We have constructed steady-state models for white dwarfs 
accreting matter of solar composition. The steady-state model 
consists of a C-O core surrounded by a hydrogen-rich enve- 
lope of the solar abundances. At the bottom of the envelope, 
hydrogen is burned at the same rate as the star accretes it; i.e.. 



Ln=XQM, 



(1) 



where L„ is the luminosity due to hydrogen burning, X 
is the hydrogen mass fraction in the accreted matter, Q = 
6.4 X 10"* erg g"' energy generated when 1 gram of hydro- 
gen is converted to helium, and M is the mass accretion rate. 
The nu clear reaction rates are taken from'Caughl an & Fowler! 
(1 19881) with the weak and intermediate screening factors by 
iGraboske et al.l (1913). The helium layer and helium shell 
burning are neglected for simplicity, because the stability of 
hydrogen burning in a thin shell is mainly determined by the 



struct ure of hydrogen-rich layer ab ove the hydrogen burning- 
shell (ISugimoto & Fujimotolll978h . Also the helium layer is 
thin and the energy release is only 10 % of hydrogen burning 
so that the effects on the stability of hydrogen burning shell 
are negligible (see P aczyriski I (1 1 983) for the study of the ef- 
fect of heat flux from the core). The compressi onal heating 
due to accretion is included in the same way as in lKawai et al.l 
(Il987l) ; i.e., the heating rate per unit mass, Cg, is given by 



M ds 
-T- 



M dlnq' 



(2) 



where s is the entropy per unit mass, and q = Mr/M with M,- 
being the mass included in the sphere of radius r. The chem- 
ical composition of the C-O core is adopted as {Xq^Xq^Xz) = 
(0.48,0.50,0.02) where denotes the mass fraction of heavy 
elements scaled by the solar composition. Opacity is obtained 
from OPAL opacity tables (Iglesias & Rogers 1996). 

Tables 1-4 summarize the properties of several steady state 
models for given M and M: the mass of the H-rich envelope 
AMenv, '"H, 7h, Pu, Ph, and £„ at the bottom of the H-rich 
envelope, and Ln, L, R, and Teff at the surface. 

3. THERMAL STABILITY OF MODELS 

For each steady-state model, we have examined the ther- 
mal stability against a linear perturbation. Assuming that per- 
turbations to stellar structure occur without disturbing hydro- 
static balance and that they are spherically symmetric, we 
have examined thermal stability by slightly modifying the 
Henyey-type relaxation code that we use to compute stellar 
structure. As well known, in the Henyey relaxation, a correc- 
tion 5yi (i= 1,2,3,4) to a variable y, is obtained by solving 
inhomogeneous linear equations expressed as 

HSy = D (3) 
where 7^ is a 4N x 4N matrix with being the number of 
mesh points, Sy a vector consists of Syi at each mesh point, 
and D is a vector with 4N components representing the devia- 
tion of the variables y, from the correct values. In thermal sta- 
bility analysis for a stellar model, we regard 6yi as a perturba- 
tion to the variable y, which satisfies the differential equations 
for stellar structure so that the inhomogeneous terms dimin- 
ish; i.e., D = Q. We express the temporal variation of perturbed 
quantities by exp(o'f ) so that we replace the time derivative of 
entropy dSs/dt with aSs. 

We write the modified Henyey matrix as Ti. Then, we ob- 
tain homogeneous equations 

my = 0, (4) 
which govern the perturbation and an eigenvalue a. 

The eigenvalue a is obtained by searching for a value which 
makes the matrix Ti, singular If there is a positive ct, the struc- 
ture is thermally unstable with a growth rate of Tg = l/cr. If 
all the eigenvalues are negative, on the other hand, we judge 
that the stellar model is thermally stable. The obtained Tg and 
the stability for several models are summarized in Tables 1 - 
4. 

Figure [T] shows the steady-state models in the HR diagram 
computed for C-O white dwarfs of M = 0.5 - 1 .38M0 which 
accrete hydrogen-rich matter at various rates M. The lowest 
luminosity in Figure[T]corresponds to M = IQt'^^Mq yr"'. 

As M gets larger along a sequence of models with given 
mass M, the luminosity asymptotically approaches a limit- 
ing value of Lrg, which corresponds to the red-giant lumi- 
nosity determined by the "core mass - luminosity" relation 
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TABLE 1 

Steady Hydrogen Burning Models (M = 1 35Mq ) 



M(M0yr') 


1.6E-9 


1.6E-8 


1.6E-7 .. 


3.5E-7 


6.0E-7 .. 


*1.6E-7'' 


AMe„v/M0 


3.7E-7 


2.2E-7 


1.4E-7 


.. 1.4E-7 


2.2E-7 


.. *1.0E-5 


log rn/RQ 


-2.488 


-2.476 


-2.463 


.. -2.453 


-2.441 


.. *-2.456 


log rH(K) 


7.67 


7.82 


7.98 


8.05 


8.11 


*8.41 


log pH(g cm-^ ) 


2.78 


2.36 


1.83 


1.56 


1.17 


*3.22 


log PH(dyn cm"^) 


18.59 


18.31 


18.04 


.. 17.98 


17.98 


*19.83 


log en(ergs g"'s"') 


9.32 


10.54 


11.72 


.. 12.08 


12.23 


*9.34 


log L„/L0 


2.048 


3.049 


4.049 


.. 4.389 


4.621 


*4.049 


log L/Lq 


2.072 


3.073 


4.077 


.. 4.423 


4.690 


*4.087 


log R/Rq 


-2.472 


-2.456 


-2.421 


.. -2.376 


-1.746 


.. *-2.334 


log reff(K) 


5.516 


5.758 


5.992 


.. 6.056 


5.807 


*5.950 


log rg(s) 


6.86 


5.94 


5.64 


stable 


stable 





" The asterisk indicates the quantities for the model whose AMenv is artificially set to be 
10"^ Mq to mimic the "surface hydrogen burning" model. 



TABLE 2 

Steady Hydrogen Burning Models {M = 1 .25Mq ) 



M(M0 yr"') 


l.OE-9 


l.OE-8 


l.OE-7 


2.1E-7 .. 


4.5E-7 


5.35E-7 


AMenv/M0 


1.6E-6 


9.3E-7 


5.7E-7 


5.4E-7 


.. 6.8E-7 


l.lE-6 


log m/RQ 


-2.306 


-2.298 


-2.282 


-2.273 


.. -2.257 


-2.250 


log rH(K) 


7.61 


7.75 


7.91 


7.97 


8.03 


8.06 


logpH(g cm"3) 


2.72 


2.30 


1.79 


1.58 


1.20 


1.04 


log PH(dyn cm"^) 


18.46 


18.19 


17.89 


17.82 


.. 17.76 


17.78 


log en (ergs g"'s"') 


8.51 


9.74 


10.93 


11.30 


11.54 


11.61 


log L„/L0 


1.845 


2.845 


3.845 


4.167 


.. 4.499 


4.574 


log L/Lq 


1.868 


2.869 


3.872 


4.196 


.. 4.538 


4.628 


log R/Rq 


-2.283 


-2.267 


-2.229 


-2.194 


.. -2.025 


-1.511 


log 7;ff(K) 


5.370 


5.612 


5.844 


5.908 


.. 5.909 


5.674 


log Tg(s) 


7.62 


6.68 


6.14 


7.36 


stable 


stable 



TABLE 3 

Steady Hydrogen Burning Models (M = I.OM0) 



M(M0yr'') 


l.OE-9 


l.OE-8 


l.OE-7 


1.3E-7 .. 


2.5E-7 


3.6E-7 


AMenv/M0 


8.0E-6 


4.7E-6 


3.2E-6 


3.2E-6 


.. 3.9E-6 


6.7E-6 


log ni/RQ 


-2.099 


-2.087 


-2.059 


-2.053 


.. -2.034 


-2.018 


log 7k(K) 


7.57 


7.71 


7.86 


7.88 


7.93 


7.96 


log pH(g cm"') 


2.51 


2.08 


1.54 


1.46 


1.20 


0.97 


log Pii(dyn cm"^) 


18.22 


17.93 


17.61 


17.58 


.. 17.50 


17.47 


log en(ergs g"'s"') 


7.82 


9.04 


10.20 


10.32 


.. 10.59 


10.66 


log Ln/LQ 


1.845 


2.845 


3.845 


3.959 


.. 4.243 


4.401 


log L/Lq 


1.868 


2.869 


3.871 


3.986 


.. 4.274 


4.439 


log R/Rq 


-2.055 


-2.026 


-1.946 


-1.923 


.. -1.792 


-1.358 


log Tef,(K) 


5.256 


5.492 


5.702 


5.720 


.. 5.726 


5.550 


log rg(s) 


8.31 


7.39 


7.20 


8.08 


stable 


stable 



mated by 

Mstabie = 3.066 X 10"^(M/MQ-0.5357)MQyr-', (5) 

and 

Mrg = 6.682 X 10~\M/Mq - O.4453)M0yr-' . (6) 

The ratio between these two rates only slightly depends on M 
asMRG/Mstabie = 2.3,2.5,2.7, and 3.0 for M= 1.38, 1.25, 1.0, 
and 0.7 Mq, respectively. The loci on the HR diagram and 
the stability of the steady-state models agree with the results 
of Sienkiewicz (1975, 19 8^. The above r esults have recently 
been confirmed by lShen & Bildstenl ('2007). who have applied 
the analytic approach with the one zone model and examined 
the metallicity dependence. 



iPaczviiskil (Il970h . (Here the core mass corresponds to the 
white dwarf mass M.) The exact value of Lrg as a function of 
M depends on the opacity and abundance of the accreted mat- 
ter, so that we give Lrg in Table 5 from our asymptotically 
obtained values. 

As the luminosity approaches Lrg the radius of the white 
dwarf gets larger. The open and filled circles indicate ther- 
mally unstable and stable models, respectively. Models 
brighter than the critical luminosity, Lstabie, around the "knee" 
of the model sequence are thermally stable. 

Table 5 summarizes the critical luminosities Lstabie and 
A^stabie, and the limiting luminosities Lrg and corresponding 
Mrg. The critical Mstabie and Mrg are, respectively, approxi- 
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TABLE 4 

Steady Hydrogen Burning Models (M = O.SMq ) 



M(M0yr') 


l.OE-9 


l.OE-8 


7.5E-8 .. 


l.OE-7 


2.2E-7 


AMenv/Mg 


2.0E-5 


1.2E-5 


l.OE-5 


.. l.OE-5 


2.0E-5 


log rii/RQ 


-1.990 


-1.973 


-1.934 


.. -1.925 


-1.888 


log Th(K) 


7.56 


7.69 


7.82 


7.84 


7.89 


log pH(g cm"3) 


2.37 


1.93 


1.44 


1.34 


1.01 


log Pn(dyn cm"^) 


18.06 


17.76 


17.46 


17.42 


17.31 


log e„(ergs g"'s"') 


7.45 


8.65 


9.64 


9.76 


10.04 


log Lu/Lq 


1.845 


2.845 


3.720 


.. 3.845 


4.188 


log L/Lq 


1.867 


2.869 


3.746 


.. 3.871 


4.218 


log R/Rq 


-1.921 


-1.874 


-1.751 


.. -1.708 


-1.255 


log reff(K) 


5.189 


5.416 


5.574 


.. 5.583 


5.444 


log Tg(s) 


8.68 


7.83 


8.47 


.stable 


stable 



TABLE 5 

Stability Boundary for Steady Burning Models 



M/Mq 


Stability Boundary 


Red Giant Boundary 


Mstable(M0yr ') 


logLstable/iQ 


MRG(MQyr ') 


log Lrg/Lq 


0.5 


1.65E-8 


3.08 


6.9E-8 


3.71 


0.6 


3.25E-8 


3.38 


1.2E-7 


3.95 


0.7 


5.35E-8 


3.60 


1.6E-7 


4.08 


0.8 


7.75E-8 


3.76 


2.3E-7 


4.23 


0.9 


1.05E-7 


3.89 


2.95E-7 


4.35 


LO 


1.35E-7 


4.00 


3.7E-7 


4.45 


LI 


1.65E-7 


4.09 


4.4E-7 


4.53 


L2 


1.95E-7 


4.16 


5.1E-7 


4.60 


L25 


2.15E-7 


4.20 


5.4E-7 


4.64 


L30 


2.35E-7 


4.25 


5.8E-7 


4.67 


L35 


2.55E-7 


4.28 


6.0E-7 


4.70 


1.38 


2.75E-7 


4.31 


6.2E-7 


4.71 




Fig. 1 . — Loci of steady models in the HR diagram. Each sequence corre- 
sponds to a mass of white dwarf which accretes matter of solar composition 
at various rates. The accretion rate determines the luminosity because hydro- 
gen is assumed to bum at the same rate of accretion. Open and filled circles 
indicate thermally unstable and stable models, respectively. 

Figure |2] shows the hydrogen-rich envelope mass AMgnv as 
a function of M for each mass of the accreting white dwarf. 
The symbols have the same meanings as those in Figure [T] 
The larger the white dwarf mass, the smaller the envelope 



Fig. 2. — The mass of the hydrogen-rich envelope versus accretion rate for 
various white dwarf masses of steady models. Open and filled circles indicate 
thermally unstable and stable models, respectively. 

mass for a given accretion rate. The mass AMgnv attains min- 
imum at the accretion rate Mstabie in equation (3) correspond- 
ing to Lstabie- FoT Smaller AMgnv, 7h would be too low to give 
enough L„. ForM > Mstabie, the models are thermally stable. 
Figure |3] shows growth time (= 1 /a) of the thermal insta- 
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Fig. 3. — Growth time of the thermal instabihty of the hydrogen burning 
shell is plotted against the accretion rate for each white dwarf mass. 

bility of steady-state white dwarf models as a function of the 
accretion rate M for each white dwarf mass M. For given M, 
the growth time decreases as M increases when M is relatively 
low. As M increases, the growth time reaches a minimum and 
increases rapidly toward the stability boundary. For given M, 
the growth time is longer for smaller M. 

To interpret the above results, we follow the description 
of basic physics of the stability of nuclear burning in a 
thin shell as summarized in Sugimoto & Fujimoto (1978) and 
ISugimoto & Miyajil (Il981h . If a temperature perturbation of 
dT{Mr) > is applied to a nuclear burning shell, it leads 
to the positive change in specific entropy of (5s(M,-) > in 
most cases mainly because of the larger temperature sensitiv- 
ity of nuclear reaction rate over radiative loss. Then 6s{M,.) 
induces the hydrostatic readjustment and the resultant tem- 
perature change is given as: 



S\n T = — 5s, 

c* 



cp 



d\nT\ 6\nP 



d\nP 



6s 



(7) 



(8) 



Here cp and c* denote the ordinary thermodynamic specific 
heat and the gravothermal specific heat, respectively. 

If c* > 0(< 0), Sin r > 0(< 0) for Ss > and shell burn- 
ing is thermally unstable (stable). Thus the stability is deter- 
mined by the hydrostatic readjustment 5ln P/5s(< 0) due to 
expansion (for 6s > 0), which depends on on the following 
two factors. 

(1) Geometry: In the extremely thin shell, the pressure at the 
bottom of the thin envelope is determined as: 



GM-(1- 



■q) 



(9) 



where R is the radius of the white dwarf, and P is determined 
only by the column density above the radius r. Therefore, the 
effect of expansion for 6ln P/6s is too small to cool the shell 
and to stabilize nuclear burning. This is the main reason for 



Fig. 4. — The properties of H-buming shells in accreting white dwarfs are 
shown in the white dwarf mass M - accretion-rate M plane. If the accretion 
rate is lower than the solid line of Mstablc. H-burning shells are thermally 
unstable. A dashed line traces the locus of the envelope mass AMsm/MQ. 
For given M and M, the envelope masses of these steady-state models are 
s maller than the e nvelope masses of the "ignition" models given in Figure 9 
in lNomotol 119821) because of the higher entropy in the steady-state models 
than th 'ignition" models (see text for more details). In the area between the 
solid (Mstablc) and dash-dotted (Mrg) lines, the H-buming shell burns steadily 
and the star is located around the "knee" or the horizontal branch on a locus 
of steady-state white dwarf models. Above the dash-dotted line for Mrg, 
the stellar envelope is expanded to a red-giant size and a strong wind occurs. 
Dotted line indicate the Eddington accretion rate Mem as a function of M. 

thin shell burning to be unstable, being the case for low M. 
For high M, entropy at the burning shell is larger, thus leading 
to a more extended envelope as seen in Figure [T] Then the 
effect of hydrostatic readjustment (expansion) is larger and 
tends to stabilize nuclear shell burning. 
(2) Equation of State: If electrons are degenerate, P depends 
only weakly on T, which makes the effect of expansion too 
small to stabilize nuclear burning upon 6T > 0. This is the 
case for low M and thus low s at the burning shell. On the 
contrary, for high M and L, especially, near their RG values, 
radiation pressure is important and its large T s ensitivity sta- 
bilizes shell burning (see also dShen & Bildstenll2007i) ). 

At Lstabie <L <Lf.G (or Mstabie < M < Mrg), therefore, 
these combined effects of radiation pressure and the extended 
envelope structure lead to stable burning. This is why MstaWe 
is smaller than M rg by only a factor of ^ 2.3 - 2.7. (See 
'Fuiimoto'(1982b) for further details on the stability of the thin 
shell approximation models.) 

Figure m summarizes the properties of steady-state models 
accreting hydrogen-rich matter. The vertical axis is the accre- 
tion rate M and the horizontal axis is the white dwarf mass M. 
The model properties are classified as follows: 
(1) The models are thermally unstable in the area below the 
solid line to show Mstabie- The dashed line indicates the loci 
where the envelope mass AMenv is constant; the value at- 
tached to each dashed line indicates AMenv/A^Q. These en- 
velope masses of the s teady state rnodels tend to be smaller 
than those obtained bv lSienkiewiczl d 1975b by a factor of 1.2 
- 1.3, but can be regarded as being consistent in view of the 
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updated opacity in the present study. 

(2) In the area above the solid hne of Mstabie, accreting white 
dwarfs are thermally stable so that hydrogen burns steadily in 
the burning shell. 

(3) Above the dash-dotted line forMRc, the accreted matter is 
accumulated faster than consumed into He by H-shell burn- 
ing. As a result the a ccreted matter is pile d up to form a 
red-giant size envelope ' Nomoto et al.l (Il979l) . The accretion 
rate is limited by the Eddington's critical accretion rate MEdd 
which is defined as 

47rcR 

Mem= . (10) 

K 

For K = 0.2(1 +X), Mem is shown by the dotted line in Fig- 
ure |4] Note that, at this critical rate, the luminosity produced 
by the gravitational energy release associated with accretion 
reaches the Eddington's critical luminosity^ 

This diagram is similar to Figure 9 in lNomotol(ll982h ex- 
cept tha t the envelope m ass in the latter is for the "ignition" 
models . ^N omotol (1 1 982h obtained the "ignition" model by cal- 
culating the time-dependent evolution of the mass accreting 
white dwarf with the Henyey code to find the stage where 
the nuclear energy generation rate e,, first exceeds the ra- 
diative energy loss rate a t the bottom of the envelope (see 
dTownslev & Bildstenll2004i) for the most updated "ignition" 
mass in one zone models). The envelope mass of the "igni- 
tion" model thus determined is larger than the envelope mass 
AMenv of our steady-state models because of the lower en- 
tropy in the "ignition" model for given M and M. 

The stability of our steady-state models are consistent 
with the previous computations for long-t erm evolu ti ons of 
accreting white dwarfs. For example, ISion et al.l (Il979h 
have found that a 1.2M0 white dwarf accreting at a rate 
1.03 X 10~^Mq yr"' gives rise to repetitive hydrogen shell 
flashes, while a 1.3M0 white dwarf accreting at a rate 
2.71 X lO'^M^ yr~' undergoes stable hydrogen burning. 
iPaczviiski & Zvtkowl (11978') have also shown that, for a 
O.8M0 white dwarf, the stability boundary of the hydrogen 
burning shell is located around M 10~^Mq yr~'. Consulting 
with Figures |2] and m we can confirm that those evolutionary 
results agree very well with our results for steady-state mod- 
els. 

4. DISCUSSION 
4.1. On "Surface Hydrogen Burning" Models 

IStarrfield et al.l (l2004t) wrote that their "surface hydrogen 
burning" models of mass accreting 1.25M0 and 1.35M0 
white dwarfs are thermally stable for accretion rates ranging 
from 1.6 X IO-'^Mq yr"' to 8.0 x lO'^M© yr"'. The stabiHty 
property of their models, however, differs from those of our 
steady-state models as wel l as all other models (Sienkiewicz 
[19801: iFuiimotol Il982ai: IShen & BildstenI I2007h . Our mod- 
els indicate that the hydrogen burning shell in the 1.35M0 
model is thermally unstable if the accretion rate is less than 
2. 5 x 10"^M(n yr" ' (2.1 X IO'^Mq yr"' for 1.25Mq) (Fig.©. 

IStarrfield et al.l (l2004h also wrote that the mass accretion 
onto the hot white dwarf just after a nova explosion leads 
to a stable surface hydrogen burning. However, the time- 
dependent calculations by Pri alnik & Kovetzl (Il995h indicated 
that the H-accretion onto the hot white dwarfs with the inte- 
rior temperatures of 1 - 5 x 10^ K leads to a shell flash for 
M = 0.65- 1.4 Mq. 

In addition to the difference in the stability property, radii 
of Starrfield et al.'s models tend to be smaller than those of 



our models (Fig. |5] below). In particular, their highest lumi- 
nosity model (M = 1.35M0 and M = 8 x 10~''Mq yr"') has a 
white-dwarf size, while our results indicate that such a high 
accretion rate makes the star of a red-giant size (Fig.|4]l. 

These discrepancies are caused by the extremely 
coarse zoning adopt e d in Starrfield et al.'s computa- 
tions. IStarrfield et al.l (120041) adopted a surface zone mass 
of 1O~^M0, which is much larger than the entire envelope 
mass of the steady-state models of M = 1 .35M0 and 1 .25Mq 
as seen in Figures |2] and |4] This means that the envelope 
of the "surface hydrogen burning" model is approximated 
by a single zone having a single temperature and density. 
Furthermore, the "surface zone" is much deeper than the 
realistic hydrogen-rich envelope of the steady-state model 
corresponding to the same white dwarf mass and the accretion 
rate. 

Table 1 compares the two white dwarf models with M = 
I.35M0 accreting hydrogen-rich matter at a rate of M = 
1.6 X IO'^Mq yr"'. The steady-state model calculated in the 
present study has AMenv = 1 .4 x 10"^Mq and log Tn(K) = 7.98 
(4th column). For the model with the asterisk (7th column), 
the mass of the hydrogen-rich envelope is artificially set to 
be AMenv = IQ~^M(7), which is the same as the "surface zone 
mass" adopted by Starrfield et al. (2004). According to equa- 
tion (|9]l and Tjj oc f h AMenv, the temperature at the burn- 
ing shell (log7H(K) = 8.41) is much higher than those of the 
steady-state model. 

Such a high temperature is comparable to that of the sur- 
face zone of Starrfield et al. (2004), from which we see the 
reason why iStarrfield et al.l (120041) obtain very high temper- 
ature at the "surface zone". They treated the envelope be- 
tween the region of log(l -q) (-5)-(-22) by a single mass 
zone, while our steady-state models resolve the H-rich enve- 
lo pe with ^ 50 mass zon es. Obviously, the zoning adopted 
by St arrfield et al.l (120041) is too coarse to obtain a physically 
realistic stellar model. 

In the heavy envelope model, the temperature at the hy- 
drogen burning shell is so high that all accreted hydrogen 
burns in one typical t ime step to compute mass accretion, as 
IStarrfield et al.l (|2004|) states "it takes less time than the time 
step (~ 2 X 10^ s) for all the infalling hydrogen to burn to 
helium in this zone". In this case, the nuclear energy gener- 
ation rate £„ is determined not by the temperature-dependent 
nuclear reaction rate but by the supplying rate of nuclear fuel 
as 

XQM 

en= ... ■ (11) 

AMenv 

Despite such a high temperature as logr(K) = 8.41, the en- 
ergy generation rate thus determined is Cn = 2.2 x 10^ ergs 
g"' s"', which is much lower than the /3-limited reaction rate 
of the hot CNO cycle e^j = 6 x 10'^(Xcno/0.01) ergs g"' s"'. 
Because XQM / AMf^nv is constant, being independent of the 
temperature, the nuclear burning is stable; it is also steady as 
expressed by equation (1). In other words, the assumed enve- 
lope mass AMenv is too large and hence the temperature at the 
nuclear burning shell is too high for the mass accretion rates 
they assumed. All the accreted hydrogen-rich matter should 
have been consumed long before it is pushed into a layer as 
deep as M-M, - IO'^Mq (iNariai et al.lll980h . 

4.2. Comparisons with Supersoft X-ray Sources 

Supersoft X-ray sources are suggested to be accreting white 
dwarfs in which steady hydrogen burning is taking place (e.g.. 




log T^ff (K) 

Fig. 5. — Loci of steady-state models in the HR diagram togetlier with Starrfield et al.'s (2004) results and positions of several supersoft X-ray sources. The 
filled small circles on the solid line are the same solutions as in Figure [T] The white dwarf mass is attached to each curve. The solid line-parts indicate the stable 
steady-state solution, while the dotted line parts correspond to the unstable solutions as shown in Figure[T] In the dashed line parts, the optically thick winds are 
predicted to blow ( Kato & Hachisu 1994), so that supersoft X-rays could be absorbed by the wind matter and not be detected. The solid line with the crosses show 
the positions of the "surface hydrogen burning" models (1.35 Mq) (Starrfiel d et al.i,2004i) . Positions of several supersoft X-ray sources are plotted by squares 
(taken from Figure 1 of Stairfield et al. 2004), i.e., R09 (RX J0925.7-4758), C87 (CAL87), C83 (CAL83), IE (1E0035.4-7230), ROO (RX J0019.84-2156), R04 
(RX J0439.8-6809). In addition, the position of the supersoft X-ray source RX J0513.9-6951 (R05) is also plotted. 



Ivan den Heuvel et all Il992h . IStarrfield et al] (120041) argued 
that the properties of their "surface hydrogen burning" models 
agree with some of supersoft X-ray sources. In this subsec- 
tion, therefore, we compare our steady burning models and 
the "surface hydrogen burning" models in the HR diagram 
with several observed supersoft X-ray sources. 

In the HR diagram of Figure |5] we show the loci of the 
steady burning models for given M. The filled small circles 
on the solid lines indicate the thermally stable models in Fig- 
ure [U while the dotted lines correspond to the unstable solu- 
tions. For comparison , the "surface hydrog en burning" mod- 
els with M = 1 .35Mq (IStarrfield et al.ll2004t) are shown by the 
crosses. In contrast to the "surface hydrogen burning" models, 
the sequences of stable models shown by the solid lines bends 
rightward because of the increasing radius with increasing M. 

The large boxes in Figure[5]indicate the positions of several 
supersoft X-ray sources taken from IStarrfield et alj (l2004t) . 
The bolometric luminosities of supersoft X-ray sources are 
very uncertain because very soft X-rays are easily absorbed 
by neutral hydrogen. For example, although the luminosity 
of CAL 87 in LMC is lower than the 1.38Mq white dwarf 
sequence in Figure |5] the bolometric luminosity of CAL 87 
has been esti mated to be a fe w times lO^'^ergs s ~', far above 
the error box (iGreine r 2000; G reiner et al.ll2004l) . Taking into 
account such large uncertainties in the bolometric luminos- 
ity, we may consider that the loci of luminous supersoft X-ray 
sources in the HR diagram are not inconsistent with the loci 
of thermally stable models (solid line parts). We suggest that 



no supersoft X-ray source is expected in the dashed line re- 
gion, becau se optically thick winds are predi cted to blow in 
this region (iKato & Hachisulll994tlKatoll 19961) and the matter 
in the winds would absorb most of the supersoft X-rays. 

We have added one more supersoft X-ray source in the 
LMC RX JOS 13.9-6951. ISchaeidt. Hasinger. & Trumpet 
(1 19931) estimated the surface temperature Tph ~ 30-40 eV 
and the total luminosity Lph ^ Lx ~ 2 x 10^^ ergs s~' by a 
black-body fitting. On the other hand, iGansicke et al.l (1 19981) 
obtained Lph Lx (3 - 9) x 10"^^ ergs s"' by their model at- 
mosphere of the white dwarf. This is a quasi-regular transient 
X-ray source with a relatively short X-ray-on 40 days) and 
along X-ray-off 100 days) states. iHachisu & Katol ( l2003al) 
proposed a transition model, in which the white dwarf ex- 
pands or contracts intermittently. In this model, an optically 
thick wind occurs only when the star is expanded, and the X- 
ray-off state is interpreted as this expanded phase where the 
supersoft X-rays are absorbed by the wind materials. 

The existence of such a transient supersoft X-ray source as 
shown in Figure |5] may support our results that there exist en- 
velope expansions for high mass accretion rates. Note that in 
Starrfield et al.'s models the envelopes do not expand even for 
higher mass accretion rates. This transient X-ray phenomenon 
is common at least in the LMC and our Galaxy, because such 
a system exists not only in the LMC (in a lower metallicity 
environment) but also in our Galaxy. A Galactic supersoft X- 
ray source, V Sge, which is not sho wn in Figure [5] is a sister 
system to RX J 0513.9-6951 (e.g.. IGreiner & van TeeselingI 
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fT998tlHachisu & Katoll20031j) . 

The above comparison shows that the properties of lumi- 
nous supersoft X-ray sources are consistent with our steady 
hydrogen-burning white dwarf models. 

5. CONCLUDING REMARKS 

We have reinvestigated the properties of accreting white 
dwarfs by constructing steady-state models, in which hydro- 
gen shell burning consumes hydrogen at the same rate as the 
white dwarf accretes it. We have confirmed that these steady 
models are stable only when the accretion rate is higher than 
^stable in equation (3). 

Our results contradict the "surface hydrogen burning" mod- 
els by Starrfield et al. (2004) who found quiescent stable hy- 
drogen burning for a wide range of accretion rates down to 
M - IO-'^Mq yr-i . We believe we have shown that quiescence 
of "surface hydrogen burning" results from the too large zone 
mass (~ IO'^'Mq) in the outer part of the numerical models, 
and that hydrogen burning must occur in a much more super- 
ficial layer (- lO'^M©). 



We have shown that the positions on the HR diagram of 
most of the luminous supersoft X-ray sources are consistent 
with the white dwarfs accreting matter at rates high enough 
for hydrogen shell burning to be thermally stable. 

The narrow range of the accretion rate that leads to 
the steady and stable hydrogen burning has been con- 
firmed. Despite that property, how the progenitor white 
dwarfs for SNe la could grow their masses to the Chan- 
drasekhar m ass by accretin g hy droge n -rich matter has been 

(e.g.,ll996); Li & van den Heuvell 
("e.g.. [l999a,b ) : Langer et al. (e.g., 

(e.g.. 



discussed in Hachisu et al 
1997); Hachis u et al 



I260O); Han & Podsiadlowskii re.g.. l206l : lNomoto et al 
l200a ,20051) 
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